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I .  INTRODUCTION 


An  essential  requirement  for  many  Navy  communication  and  radar 
receiving  systems  is  the  ability  to  operate  in  the  face  of  strong  inten¬ 
tional  interference  signals.  One  technique  that  has  been  shown  to  be 
capable  of  providing  a  margin  of  30+  dB  against  jamming  is  the  use  of  an 
adaptive  antenna  array.  Adaptive  arrays  form  their  beam  in  response  to 
the  signals  impinging  on  the  array;  for  defeating  jammers,  the  array 
places  pattern  nulls  in  the  directions  of  the  undesired  jamming  signals. 

Most  adaptive  arrays  resemble  conventional  deterministically 
steered  arrays  in  many  respects.  Each  element  has  a  complex  weight,  and 
the  weighted  element  outputs  are  summed  in  an  electronic  summing  element 
to  obtain  the  array  output.  In  the  case  of  an  adaptive  array,  however, 
the  output  is  processed  to  derive  control  voltages  that  are  fed  back  to 
the  complex  weights. 

In  this  report  we  describe  a  new  type  of  very  compact  adaptive 
antenna  array  that  does  not  use  an  electronic  summing  element.  This 
array,  which  we  term  a  reactively  steered  adaptive  array  (RESAA) , 
differs  from  other  adaptive  arrays  fundamentally  in  that  only  one  ele¬ 
ment  is  physically  connected  to  a  receiver.  The  remaining  elements  are 
parasitic,  and  the  pattern  is  formed  according  to  the  values  of  the 
reactive  terminations  on  these  parasitic  elements.  The  array  "summing" 
occurs  in  the  currents  on  the  one  element  connected  to  the  receiver; 
thus,  the  electronic  summing  device  is  eliminated. 

This  report  discusses  experimental  results  and  the  supporting 
theory  for  a  RESAA  listing  of  five  microstrip  rectangular  patch 
elements.  The  array  operates  at  a  frequency  of  4,0  GHz.  This  array  is 
intended  primarily  as  a  generic  form  of  a  RESAA  to  demonstrate  the 
capabilities  of  the  technique  and  to  use  as  a  test  bed  for  development 
of  an  automatic  control  method.  However,  because  of  their  low  profile, 
microstrip  elements  are  increasingly  being  favored  for  flush-mounted 
conformal  antennas  on  missiles  and  aircraft.  The  application  of  adap¬ 
tive  array  techniques  to  microstrip  elements  is  of  particular  interest, 
because  the  mitigation  of  jammers  using  adaptive  beam-terming  tech¬ 
niques  helps  to  overcome  the  vulnerability  to  jamming  that  the  rela¬ 
tively  narrow  bandwidth  of  microstrip  elements  permits. 

This  report  is  organized  as  follows.  The  next  section  describes 
the  experimental  results  obtained  using  manual  adjustment  of  the 
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reactive  terminations  on  a  five-element  microstrip  array.  Section  III 
presents  the  theory  describing  the  operation  of  a  RESAA.  This  theory 
is  used  to  interpret  the  experimental  results  of  the  microstrip  arrav , 
but  a  detailed  comparison  of  theory  and  experiment  is  not  attempted. 
Rather,  the  theory  is  used  in  Section  IV  to  analyze  one  possible  auto¬ 
matic  adaptive  control  method  that  could  be  used  with  any  RESAA.  In 
Section  V  a  discussion  is  given  of  future  plans  and  recommendations. 


II.  EXPERIMENTAL  MICROSTRIP  ARRAY  RESULTS 


A.  ARRAY  DESIGN 

To  demonstrate  the  RESAA  concept,  the  microstrip  array  shown  in 
Figure  1  was  fabricated  and  tested  using  manual  control  of  the  phase 
shifter  reactive  loads.  The  elements  were  designed  to  have  a  center 
frequency  at  4.0  GHz,  a  frequency  chosen  as  a  reasonable  tradeoff 
between  a  convenient  antenna  size,  fabrication  ease  of  phase  shifters 
and  antenna  elements,  and  availability  of  an  anechoie  chamber  for  test¬ 
ing.  The  spacing  of  0.75  cm  corresponds  to  a  separation  of  0.1\  at 
4.0  GHz;  the  experimental  coupling  results  of  dedlicka  er  a/.,1  along 
with  extrapolation  of  our  earlier  experience2  with  arrays  at  20  MHz, 
indicated  that  this  separation  would  produce  a  sufficient  amount  of 
mutual  coupling  between  elements.  The  substrate  material  is  Rexolite 
1422  with  a  thickness  of  1.58  mm. 


B.  PHASE  SHIFTER  AND  REACTIVE  LOAD  DESIGN 

A  variable  reflection  phase  shifter  terminating  a  transmission 
line  with  a  characteristic  impedance  Z  produces  a  variable  reactance  X 
according  to  the  equation 


X  =  Z  tan<f> 
o 


(1) 


where  (j>  is  the  phase  shift.  Microstrip  phase  shifters  were  designed 
and  fabricated  to  provide  the  necessary  amount  of  phase  shift,  using  a 


1  R.  P.  Jedlicka,  M.  T.  Poe,  and  K.  R.  Carver.  "Measured  Mutual 
Coupling  Between  Microstrip  Antennas,"  IEEE  Trans.  Antennas  Pr'opag .  , 
Vol .  AP-29  (Jan  1981),  pp.  147-149. 

2  Naval  Research  Laboratory.  "A  Compact  HE  Antenna  Array  Using 
React ively  Terminated  Parasitic  Elements  for  Pattern  Control,"  by 

R.  J.  Dinger  and  W.  D.  Meyers.  Washington,  DC,  NRL,  May  1982.  (NRL 
Memo  Report  4797,  publication  UNCLASSIFIED.) 
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FiCFKK  1.  llxper  imental  Kl'SAA  and 
Diagram  of  Antenna  Coordinate  System. 

1’  denotes  the  field  observation  point. 

varaetor  diode  mounted  between  a  'i()-olmi  microstrip  line  and  the  ground 
plane.  Figure  2  is  a  photograph  of  the  phase  shifter,  and  Figure  1 
displavs  the  phase  shift  versus  control  voltage  characteristic  of  three 
of  tlie  finished  devices.  The  control  characteristic  is  clearlv  not 
linear,  but  a  linear  character  1st  ic  is  not  required  for  this  arrav 
because  it  is  adaptive  (i.e.,  precise  knowledge  of  the  conversion 
between  control  voltage  and  phase  shift  is  not  needed  in  terming  the 
antenna  pattern).  However,  better  linearity  than  is  evident  in  Figure 
would  be  desirable  for  best  performance,  since  most  of  the  phase  control 
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FIGURE  2.  Microstrip  Phase  Shifter. 


FIGURE  3.  Reflection  Phase  Shift  as  a 
Function  of  Varactor  Bias  Voltage  for 
Three  Units. 
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occurs  only  in  about  6  volts  of  the  30-volt  bias  range.  A  linear 
control  characteristic  requires  at  least  two  varactor  diodes,'1  a  compli¬ 
cation  that  was  felt  unwarranted  for  this  proof-of-concept  array.  The 
reflection  loss  for  the  phase  shifters  was  typically  less  than  1.5  dB. 

Figure  4  is  a  photograph  of  the  rear  of  the  array,  showing  the  four 
phase  shifters  connected  to  the  feeds  of  the  four  parasitic  elements. 


FIGURE  4.  Rear  View  of  Experimental 
RESAA  Showing  Phase  Shifters  in 
Place. 


C.  ARRAY  PERFORMANCE 

A  block  diagram  of  the  installation  of  the  array  in  an  anechoic 
chamber  is  shown  in  Figure  5.  All  antenna  patterns  shown  in  this 
report  were  taken  in  the  H-plane  of  the  array.  The  measurements  pro¬ 
ceeded  in  three  steps.  (1)  The  array  was  rotated  in  the  0-plane  to 
place  the  incident  angle  of  the  anechoic  test  signal  in  a  direction  to 
represent  an  interference  signal.  (2)  The  array  output  (as  produced  by 
the  mixer-receiver  combination)  was  then  monitored  and  the  control 
voltages  for  the  phase  shifters  iteratively  adjusted  to  minimize  its 
value.  Typically,  about  three  passes  through  the  four  phase  shifter 
voltage  adjustments  were  enough  to  decrease  the  signal  level  to  its 
lowest  value.  (3)  The  phase  shifter  voltage  levels  were  left  untouched 
(i.e.,  the  reactive  terminations  were  "frozen"),  and  the  anechoic 
chamber  test  signal  was  then  used  to  record  the  antenna  pattern. 


3  R.  V.  Garver.  "360  Degree  Varactor  Linear  Phase  Modulator," 
IEEE  Trans.  Microwave  Theory  Tech.,  Vol .  MTT-17  (March  1969),  pp.  137- 
147. 
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FIGURE  5.  Diagram  of  the  Installation  of  the  Experimental  RESAA  in 
an  Anechoic  Chamber. 


The  results  of  this  process  are  shown  in  Figure  6.  Recorded 
patterns  are  given  for  incident  signal  angles  of  0,  30,  45,  and  60 
degrees.  The  adjustment  of  the  reactive  terminations  is  observed  in 
those  patterns  to  place  a  minimum  in  the  pattern  towards  the  direction 
of  the  interference  source.  This  minimum  has  a  width  of  about  25 
degrees  and  a  depth  (as  measured  relative  to  the  average  pattern  ampli 
tude)  of  about  30  dB.  The  other  relatively  minor  features  in  the 
patterns  result  from  unused  degrees  of  freedom;  that  is,  more  than  one 
set  of  reactive  termination  values  will  produce  the  desired  minimum  in 
the  pattern,  and  control  over  features  other  than  the  minimum  requires 
additional  information  to  choose  the  correct  set  of  termination  values 
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FIGURE  6.  Log-Polar  Plots  of  the  Measured  Radiation  Pattern 
After  a  Signal  Incident  in  the  Direction  Indicated  by  the 
Arrow  was  Minimized  by  Reactive  Load  Adjustment:  0-  and 
30-Degree  Incidence.  0  dB  is  an  incident  power  of  7.0  dBi , 
as  measured  by  a  standard  gain  horn. 
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Figaro  7  shows  the  sensitivity  of  the  ittem  minimum  to  frequency. 
In  this  figure,  an  incident  signal  at  3.95  GHz  was  nulled  and  the  reac¬ 
tive  terminations  frozen.  Patterns  were  then  measured  at  the  frequen¬ 
cies  on  both  sides  of  3.95  GHz.  Figure  7  indicates  that  as  the  frequency 
is  changed  from  that  at  which  the  null  was  formed,  the  minimum  shifts  in 
direction  and  the  depth  of  the  minimum  decreases. 


0° 


FIGURE  7.  Log-Polar  Plot  of  the  Measured  Radiation  Pattern  for  Three 
Frequencies.  The  reactive  loads  were  adjusted  at  a  frequency  of 
3.95  GHz  and  left  untouched  for  the  other  two  frequencies. 


Figure  8  is  a  plot  of  pattern  value  at  the  minimum  point  versus 
frequency;  defining  the  nulling  bandwidth  as  the  frequency  range  over 
which  tiie  null  is  at  least  20  dB  below  the  pattern  maximum,  a  value  of 
about  40  MHz  is  obtained  for  this  array. 

To  summarize  the  experimental  results,  the  following  points  are 
emphasized . 

•  The  interference  nulling  procedure  consistently  resulted  in  an 
antenna  pattern  with  a  null  that  was  narrow  (25  degrees)  and 
reasonably  deep  (at  least  30  dB) . 

•  The  null  could  be  steered  to  essentially  any  angle  of  incidence, 
except  for  grazing  incidence  on  the  antenna  ground  plane. 
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In  the  next  section,  equations  are  derived  for  the  radiation 
pattern  of  a  RESAA,  both  for  a  general  form  of  array  element  and,  speci¬ 
fically,  for  a  microstrip  patch  element;  the  pattern  is  expressed  as  a 
function  of  the  reactive  termination  loads  on  the  parasitic  elements. 

We  do  not  attempt  a  point-by-point  fit  between  the  experimentally  meas¬ 
ured  patterns  and  the  patterns  predicted  by  the  theory.  We  will  note 
that  the  general  behavior  predicted  by  the  theory  is  observed  in  the 
experimental  measurements  and  then  proceed  to  use  the  theory  to  devise 
and  evaluate  an  automatic  adaptive  control  technique  for  the  RESAA. 

III.  REACT! V ELY  STEERED  ARRAY  THEORY 


Harrington4  first  presented  the  tiieory  for  an  array  of  react ively 
steered  (wire)  dipole  antennas  and  later5  extended  this  theory  to  an 
array  of  dielectrically  loaded  waveguide-fed  slot  antennas.  All  of  his 
analyses  treated  the  array  deterministically  in  the  sense  that  reactive 
load  selection  was  accomplished  by  an  open-loop  technique.  In  this 
section,  we  extend  Harrington's  analysis  to  the  microstrip  array. 


A.  ARRAY  PATTERN  FOR  A  GENERAL  PARASITIC  ARRAY 

In  analyzing  the  wire  dipole  array,4  Harrington  used  the  Thevenin 
equivalent  circuit  of  Figure  9,  which  shows  a  load  network  and  an  N-port 
antenna  system  connected  in  parallel  with  series-driving  voltages.  The 
terminal  equation  for  the  load  and  antenna  system  of  Figure  9  is 


f°c  -  (ZA  +  zji 


(2) 


0  C  -> 

where  V  and  I  are  the  column  vectors  of  the  Thevenin  equivalent  volt¬ 
ages  and  port  currents,  respectively.  The  matrices  [ZA]  and  [ZjJ  are 
the  open-circuit  impedance  matrices  of  the  antenna  system  and  load 
system,  respectively.  The  electric  field  radiated  by  the  antenna  system 
can  then  be  written  by  superposition  as 


I  E°C 
n  n 


(3) 


4  R.  F.  Harrington.  "Reactively  Controlled  Directive  Arrays," 
IEEE  Trans.  Antennas  Pro pap.  ,  Vol.  AP-26  (May  1978),  pp.  390-347. 

5  J.  Luzwick  and  R.  F.  Harrington.  "A  Reactively  Loaded  Aperture 
Antenna  Array,"  IEEE  Trans.  Antennas  Propaa .,  Vol.  AP-26  (July  1978), 
pp.  870-873. 
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where  Ea  is  the  field  radiated  when  a  unit  current  exists  at  port  n 
and  all  otiier  ports  are  open  circuited.  Equation  2  can  be  inverted  and 
solved  for  I  and  substituted  into  Equation  3  to  give 


E  = 


EOC[Z, 


zL]-^oc 


(4) 


FIGURE  9.  Thevenin  Equivalent  Circuit  of  an  Antenna  Arrav  with 
Parasitic  Elements. 


Since  only  the  center  element  denoted  by  p  =  (N  +  l)/2  is  driven  in 
this  array,  Voc  is  a  vector  with  all  zeros  except  for  the  pc"  element. 
Hence,  Equation  4  can  be  written  as 

E  =  V  7  E  { (ZA  +  Z, }"1}  ,  (3) 

p  '  n  A  L'  np  ' 

where  (-lnp  denotes  the  np  element  of  the  matri-  \Z\  +  ZjJ  ^  and  Vp  is 
the  open-circuit  voltage  at  the  center-driven  element.  For  a  given 
type  of  antenna  element  and  array  geometry.  Equation  5  can  be  used  to 
determine  the  (unnormalized)  array  pattern  by  computing  the  elements  of 
[Z^l  and  l Zj^ ]  and  substituting  the  appropriate  form  for  E^. 
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B.  MICROSTRIP  ARRAY  PATTERN 

This  method  of  analysis  using  N-port  network  theory  clearly  is 
suitable  for  wire  antennas  and  lumped  element  loads;  it  was  used/  for 
example,  to  synthesize  antenna  patterns  at  20  MHz.  The  applicability  of 
the  theory  to  aperture  antennas  and  transmission  line  loads  is  not 
obvious;  however,  Harrington  and  Mautz6  have  shown  that  N-port  network 
theory  can  be  used  for  aperture  problems  by  dividing  the  analysis  into 
two  separate  regions  by  using  the  equivalence  principle.  One  region 
contains  only  the  loads  (i.e.,  the  space  "behind"  the  aperture),  and  the 
other  region  is  the  half-space  into  which  the  array  radiates.  Two 
aperture  impedance  matrices  are  computed,  one  for  each  region,  that  are 
independent  of  each  other.  Luzwick  and  Harrington5  used  this  approach 
to  analyze  a  reactively  steered  array  of  dielectrically  loaded  waveguide- 
fed  slot  antennas  and  arrived,  after  a  lengthy  analysis,  at  an  equation 
of  the  same  form  as  Equation  5  (actually,  the  admittance  equation  that 
is  the  dual  of  Equation  5,  since  they  chose  to  use  an  admittance  formu¬ 
lation).  Therefore,  without  proceeding  through  a  complete  derivation 
starting  from  equivalent  current  sheets  for  the  microstrip  patch  ele¬ 
ments,  we  assert  that  Equation  5  can  be  used  directly  to  write  down  the 
array  radiation  pattern.  The  task  then  is  to  derive  expressions  for  (1) 
[ZaK  Che  diagonal  matrix  of  reactive  terminations,  (2)  [Zl]»  the  mutual 
impedance  matrix  for  the  radiation  half-space,  and  (3)  En,  which  is 
equivalent  to  the  radiated  electric  field  of  an  isolated  microstrip 
patch  element. 


1 .  Load  Matrix  [Z^] 

The  load  matrix  has  the  form  (Z, ) .  .  =  Z,  .  6 .  . ,  except  that  Z,  =  0, 

L  i]  Li  ij  Lp 

corresponding  to  the  driven  element  (assumed  to  be  matched).  The  form 
of  each  load  is  given  by  Equation  1,  so  that  the  entries  along  the 
diagonal  of  [Z^]  can  be  written  compactly  as 


(Z,  ) .  .  =  jZ  tan<J> . 
L  li  o  l 


where  4>  ^ 


is  the  phase  shift  for  the 


.th 

i 


element . 


(6) 


6  R.  F.  Harrington  and  J.  R.  Mautz.  "A  Generalized  Network  Formu¬ 
lation  for  Aperture  Problems,"  IEEE  Trann.  Antannap  Propan.,  Vol.  AP-24 
(Nov  1976),  pp.  870-873. 
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2 .  Mutual  Impedance  Matrix  [2^1 

Several  techniques  for  computing  the  mutual  impedance  between 
microstrip  patch  antenna  elements  have  been  reported,7-1®  all  of  which 
require  considerable  machine  computation  and  none  of  which  can  be  writ¬ 
ten  down  in  convenient  analytical  form.  The  equations  due  to  Krowne  and 
Sindoris1®  are  the  most  straightforward;  however,  their  applicability  is 
somewhat  limited,  since  the  equations  are  apparently  valid  only  for  very 
closely  spaced  elements.  For  design  and  analysis  of  arbitrary  patch 
arrays,  theoretical  expressions  for  the  mutual  impedance  over  a  wide 
range  of  array  parameters  are  needed.  Although  we  have  done  some  calcu¬ 
lations  using  the  equations  in  Reference  10,  the  calculations  presented 
below  have  avoided  the  problem  by  using  mutual  impedance  matrix  values 
obtained  by  ’'icasu.Y'cments  on  the  experimental  antenna  shown  in  Figure  1. 
These  measurements  were  taken  using  an  automatic  network  analyzer ,  a 
discussion  of  which  is  given  in  Appendix  A. 

3.  Radiation  Field  E 

- - . — — - — — - n 

We  use  a  transmission  line  model  of  a  microstrip  patch  element7 8 * 10 11 
with  the  coordinate  system  shown  in  Figure  1.  Radiation  occurs  mainly 
from  fringing  fields  at  the  open-circuited  ends,  and  each  end  is  assumed 
to  radiate  as  a  dipole  element.  The  total  radiated  field  from  one  patch 
element  is  then  given  by  the  superposition  of  two  dipole  fields  sepa¬ 
rated  by  the  patch  width  L.  To  within  a  factor  that  is  constant  for  all 
elements,  the  radiation  pattern  for  E^  is  given  by 

2sin(k  hcos0/2)  jk  x  cosy 

(Ez)n  = - iT^0 -  cos(koLcos6/2)e  °  n  ,  (7) 

o 

where  h  is  the  substrate  thickness  and  kG  =  2n/A.  The  exponential 
factor  in  Equation  7  accounts  for  the  phase  variation  across  the  array, 
with  the  phase  center  taken  at  the  center  element. 


7  M.  Malkomes.  "Mutual  Coupling  Between  Microstrip  Patch  Antennas." 
Electron*  Lett.,  Vol.  18  (June  1982),  pp.  520-521. 

8  S.  M.  Ali,  T.  M.  Habashv,  and  J.  A.  Kong.  "Resonance  in  Two 
Coupled  Circular  Microstrip  Disk  Resonators,"  J.  J/y  /.  /’*/;< J.  ,  Vol.  53 
(Sept  1982),  pp.  6418-6429. 

8  C.  M.  Krowne.  "E-Plane  Coupling  Between  Two  Rectangular  Micro¬ 
strip  Antennas,"  ELeelpot;.  Let’-.,  Vol.  16  (Mar  1980),  pp.  211-213. 

10  C.  M.  Krowne  and  A.  R.  Sindoris.  "H-Plane  Coupling  Between 
Rectangular  Microstrip  Antennas,"  FU?e1}'o>'.  Lett.,  Vol.  16  (Mar  1980), 
pp.  211-213. 

11  R.  E.  Munson.  "Conformal  Microstrip  Antennas  and  Microstrip 
Phased  Arrays,"  I  KEF.  ''.'vane.  Antcr'hin  Pvopaij . ,  Vol.  AP-22  (  1974),  pp.  74- 
78. 
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The  array  pattern  can  then  he  written  as 


F  (0 )  = 


2sin(k  heosO/2) 
o _ 

k  licosO 
o 


eos(k  I.cosH/2)  J  l  [  2  +Z 


n=l 


,  jk  x  r os6 
,  -1 ,  on 

>  e 
np 


(8) 


By  reciprocity.  Equation  8  can  also  he  viewed  as  proportional  to  the 
array  output  voltage  when  the  array  is  receiving  a  signal  polarized  in 
the  Ez  direction  and  incident  at  an  angle  0.  The  presence  of  the  matrix 
inverse  in  Equation  8  in  general  makes  it  difficult  to  display  the 
explicit  dependence  of  F  on  the  reactive  loads.  However,  for  two-  and 
three-element  arrays,  the  inverse  can  be  computed  analytically  in  a 
reasonable  fashion;  in  Appendix  B  Equation  8  is  written  out  explicitly 
for  these  two  cases. 


The  connection  between  a  react ivelv  steered  array  and  a  conventional 
fully  driven  array  can  be  made  by  noting  that  Equation  8  can  be  written 
in  the  standard  form  of  a  fully  driven  array  as 


jk  x  cosO 

F ( 0 )  =  ge(0)  l  Ine  °  n 
n 


(9) 


in  which  g  (0)  is  the  element  pattern  and  the  feed  currents  I  are 
e  n 


I 

n 


IIZA  +  ZI. 


(10) 


By  comparison,  the  feed  currents  of  the  conventional  array  are  given 

.  ,  ,  ,conv  .  i,i  10n  ,  ,  .  ,  ,  .  ,  t  . 

simplv  bv  I  =  A  A  io  ,  where  A  is  the  complex  weight  for  the 

•  n  n  1  iv  n 

n  element.  The  most  noticeable  differences  between  the  feed  currents 

of  the  react ivelv  steered  array  and  the  conventional  array  are  the 

former's  nonlinear  relationship  between  current  and  controlled  variable 

(the  Xj/s)  and  the  dependence  of  each  of  its  currents  on  all  of  the 

controlled  variables. 

In  the  simulations  that  follow,  the  incident  signal  and  interfer¬ 
ence  are  both  assumed  to  be  Fz  components .  Additionally,  a  noncoherent 
system  with  a  square  law  detector  is  assumed,  so  that  the  array  output 
is  computed  as  V  =  [EE"]1^. 
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IV.  ADAPT  1VF.  CONTROL  SIMULATION 


A.  BACKCROUND 

The  adaptive  control  discussion  nnd  simulations  in  this  report  are 
limited  to  a  power  inversion  mode  of  operation  tor  the  array,  in  which 
the  reactive  loads  are  adjusted  to  minimize  the  incident  interference 
power.  This  mode  of  operation  was  basically  the  technique  used  in  the 
experimental  measurements  to  form  a  null. 

A  generalized  block  diagram  of  a  react ively  steered  array  connected 
for  adaptive  control  of  the  loads  is  shown  in  Figure  10.  For  a  given 


¥ 

OUTPUT 


FIGURE  10.  Diagram  of  Reactively  Steered  Array  with  a  Generalized 
Adaptive  Controller.  For  a  mode  of  operation  requiring  a  reference 
signal,  the  injection  of  the  reference  signal  is  indicated.  How¬ 
ever,  for  the  measurements  and  analvsis  in  this  paper,  no  reference 
was  used. 

angle  of  incidence  of  interference,  the  arrav  output  as  given  hv  Equa¬ 
tion  8  produces  a  surface  in  the  N-d imens ional  space  formed  bv  (N-l) 
axes  for  the  reactive  loads  and  one  axis  for  the  arrav  output.  This 
surface  is  the  error  surface  for  the  adaptive  control  algorithm;  the 
purpose  of  the  algorithm  is  to  find  the  set  of  reactive  loads  corre¬ 
sponding  to  the  (hopefullv  single)  minimum  point  of  the  error  surface. 
Most  adaptive  algorithms  (an  exception  is  random  search)  are  of  the 
steepest  descent  type  that  find  the  minimum  in  the  error  surface  by 
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beginning  at  some  point,  sousing  the  gradient  or  selected  components  of 
the  gradient  at  that  point,  and  incrementing  the  point  in  a  direction 
corresponding  to  the  steepest  value  of  the  gradient.  The  Widrov  least 
mean  squares  (LMS)  algorithm*2  in  a  full  driven  array  senses  the  error 
surface  gradient  by  correlating  the  outputs  of  each  array  element  with 
tite  total  array  output  and  then  forming  an  error  signal.  However,  for 
tlie  react ively  steered  array  only  the  output  of  the  single  center  ele¬ 
ment  of  the  array  is  available,  so  that  LMS  processing  is  not  possible. 

As  indicated  in  Figure  10,  some  type  of  modulation  of  the  reactive 
loads  must  be  provided  in  order  to  measure  the  gradient  of  the  error 
surface.  This  dither  signal  must  be  applied  and  the  resulting  influence 
on  the  RESAA  output  measured  at  the  same  time  that  the  interference  and 
desired  signals  are  being  detected. 


B.  ERROR  SURFACES 

Figures  11  through  13  present  examples  of  the  error  surface  for  a 
microstrip  RESAA.  In  these  figures  Equation  8  is  plotted  for  an  array 
with  the  following  characteristics. 

•  A  three-element  array  is  assumed  in  which  each  element  is 
identical  to  the  elements  in  the  experimental  array  shown  in 
Figure  1. 

•  Ttie  element  spacing  is  the  same  as  the  experimental  array  in 
Figure  2. 

•  The  values  of  [Z^]  in  Equation  8  are  the  values  measured  for 
the  three  center  elements  of  the  experimental  array,  as  given 
in  Appendix  A. 

For  each  error  surface  plot  the  angle  of  incidence  of  the  signal  is 
fixed,  and  Equation  8  is  plotted  as  a  function  of  X]/Zj  and  X^/Zq. 

Three  different  values  of  6  are  used  in  the  three  plots. 

The  error  surfaces  for  this  three-element  example  show  one  well- 
defined  global  minimum  that  corresponds  to  those  values  of  the  reactive 
terminations  that  reduce  the  incident  interference  to  its  lowest  value 
at  the  array  output.  The  location  of  the  minimum  changes  as  the  inci¬ 
dence  angle  changes.  In  Figure  14  the  trace  of  the  minimum  point  on  the 
X]  -  \'i  plane  is  shown. 


12  J.  E.  Hudson.  Adaptive*  Arvuit  Pvr  ru"'}  •  Ire .  London,  Institution  of 
Electrical  Engineers,  1981. 
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FIGURE  11.  Error  Surface  (Array  Output  Versus  Reactive  Load  Value) 
for  a  Throe-Element  RESAA  with  Interference  Incident  at  0  =  0  Do- 


FIGURE  12.  Error  Surface  for  a  Three-Element  RESAA  with  Inter  for 
ence  Incident  at  0  =  -45  Degrees. 
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FIGURE  13.  Error  Surface  for  a  Three-Element  RESAA  with  Interfer¬ 
ence  Incident  at  0  =  -90  Degrees. 


C.  STEEPEST  DESCENT  CONTROL  ALGORITHM 

t  h  th 

Let  X^(j)  denote  the  reactive  load  on  the  i  antenna  at  the  j 
iteration.  Then  the  discrete  steepest  descent  algorithm  states  that  the 
reactive  load  at  the  next  iteration,  X^(j+1),  is  given  by 

X.(j  +  1)  =  X(j)  -  KVxf  ,  (11) 


where  K  is  a  constant  that  controls  the  convergence  rate  and  t  is  the 
n-d imensional  error  surface.  In  terms  of  the  error  surfaces  given  in 
Figures  11  through  13,  the  reactive  loads  are  adjusted  to  follow  a  path 
that  traces  out  the  maximum  value  of  the  gradient  at  every  point  and 
that  proceeds  to  successively  lower  values  of  the  array  output. 

In  a  practical  implementation  of  steepest  descent  control,  the 
partial  derivative  in  Equation  11  can  be  measured  by  superimposing  a 
modulation  signal  AX  on  the  ractive  load  control  voltage  and  detecting 
the  corresponding  change  AR  in  receiver  output;  the  partial  derivative 
is  approximated  by  AR/AX.  Since  the  partial  derivative  is  needed  for 
each  parasitic  element,  some  technique  must  be  provided  to  separate  the 
responses  at  the  receiver  output  for  each  element.  Two  possible  methods 
for  this  separation  are  shown  in  Figures  15  and  16.  For  the  method  in 
Figure  15  the  modulation  frequency  for  each  parasitic  element  is  differ¬ 
ent,  and  synchronous  detection  is  used  to  separate  the  response  for  each 
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FI  (HIRE  1A.  Trace  of  the  Minimum  Points  in  Figures  11  Through  li  on 
a  Plane  Perpendicular  to  the  Array  Output  Axis. 


element;  this  technique  is  probably  best  implemented  using  analog  cir¬ 
cuitry.  For  the  method  in  Figure  16,  the  reactive  loads  are  modulated 
sequentially  in  time,  a  technique  that  is  probably  best  implemented 
using  digital  circuitry.  The  simulations  discussed  in  the  next  section 
applv  more  directly  to  this  latter  technique. 
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FIGURE  15.  Technique  for  Estimating  the  Error  Surface  Slope  Using 
a  Synchronous  Detector  for  Each  Element. 


FIGURE  16.  Technique  for  Estimating  the  Error  Surface  Slope  Using 
One  Synchronous  Detector  with  Time  Sharing  Between  the  Elements. 
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D.  SIMULATION  RESULTS 

The  steepest  descent  algorithm  was  simulated  on  a  Hewlett-Packard 
Model  1000  computer  assuming  the  three-element  RF.SAA  discussed  in  Sec¬ 
tion  IV-B.  To  generate  the  estimate  of  the  error  surface  slope  for 
Equation  11,  the.  change  in  the  array  output  in  response  to  a  small 
change  in  each  reactive  termination  was  determined  at  each  iteration. 
The  small  change  in  the  reactive  termination  is  referred  to  as  the  step 
size,  below. 

Figures  17  through  19  plot  the  evolution  of  the  reactive  loads  Xj 
and  X3  during  the  steepest  descent  algorithm  for  three  different  inter¬ 
ference  angles  of  incidence.  In  all  three  figures  the  initial  values 


FIGURE  17.  Variation  of  Reactive  Loads  During  Steepest  Descent 
Convergence  on  the  Error  Surface  of  Figure  11.  Interference 
incidence  angle  =  0  degrees,  step  size  =  0.01.  For  the  dashed 
curve,  K  =  0.01;  for  the  solid  curve,  K.  =  0.035. 


of  the  normalized  reactive  loads  are  taken  as  0.2,  and  curves  are  shown 
for  two  different  values  of  the  convergence  constant  K.  In  Figure  20 
the  reduction  in  the  interference  power  during  convergence  to  the  error 
surface  minimum  is  shown  for  some  of  the  cases  presented  in  Figures  17 
through  19.  Figure  21  is  a  display  of  convergence  time  (defined  as  the 
number  of  iterations  needed  to  reach  90  percent  of  the  final  value  of 
the  reactive  load)  and  the  peak-to-peak  oscillation  amplitude  of  the 
steady-state  reactive  load  value  plotted  as  a  function  of  the  value  of 
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FIGURE  18.  Variation  of  Reactive  Loads  During  Steepest  Descent 
Convergence  on  the  Error  Surface  of  Figure  12.  Interference 
incidence  angle  =  -45  degrees,  step  size  =  0.01.  For  tae  dashed 
curves,  K  =  0.01;  for  the  solid  curves,  K  =  0.028. 


FIGURE  19.  Variation  of  Reactive  Loads  During  Steepest  Descent 
Convergence  on  the  Error  Surface  of  Figure  13.  Interference 
incidence  angle  =  -90  degrees,  step  size  =  0.01.  For  the  dashed 
curves,  K  *  0.01;  for  the  solid  curves,  K  =  0.025. 
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FIGURE  20.  Change  in  Received  Interference  Signal  Power  During 
Steepest  Descent  Processing.  Solid  curve:  interference  inci¬ 
dence  angle  =  0  degrees,  K  =  0.01.  Dashed  curve:  interference 
incidence  angle  =  -45  degrees,  K  =  0.028.  Dot-dashed  curve: 
interference  incidence  angle  =  -45  degrees,  K  =  0.01. 


FIGURE  21.  Convergence  Time  and  Steady-State  Reactive  Load  .litter 
as  a  Function  of  K.  Interference  incidence  angle  =  0  degrees. 
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R.  In  Figures  22  through  24  wo  show  the  beam  history  for  several 
selected  sets  of  parameters.  Table  1  summarizes  the  parameters  used  for 
Figures  17  through  24. 


FIGURE  22.  Beam  History  During  Steepest  Descent  Processing  for  Inter¬ 
ference  Incident  at  0  Degrees. 


The  following  conclusions  can  be  drawn  from  these  results. 

•  The  interference  power  reduction  is  greatest  at  broadside  inci¬ 
dence  (30  dB)  and  is  least  at  end-fire  incidence.  The  converg¬ 
ence  time  is  related  to  the  sharpness  of  the  minimum  (i.e., 

the  second  derivatives  near  the  minimum);  the  sharper  the  mini¬ 
mum,  the  larger  are  the  values  of  the  gradient  in  Equation  11 
for  fixed  K.  Figures  11  through  13  show  that  the  minimum  is 
sharpest  at  broadside,  which  is  consistent  with  the  fast  converg¬ 
ence  . 

•  A  trade-off  between  convergence  speed  and  steady-state  fluctua¬ 
tions  in  the  reactive  load  values  can  be  made.  This  kind  of 
trade-off  is  usually  required  in  adaptive  antennas,  and  the 
reactively  steered  array  is  no  exception. 

Finally,  in  Figure  25  we  show  a  beam  history  during  steepest 
descent  processing  when  two  interference  sources  are  present,  one  at  0 
degrees  and  the  other  at  -45  degrees.  A  three-element  array  does  not 
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6  (DEGREES) 

FIGURE  23.  Beam  Historv  During  Steepest  Descent  Processing  for  Inter¬ 
ference  Incident  at  -43  Degrees. 

possess  a  sufficient  number  of  degrees  of  freedom  to  place  deep  minima 
in  the  direction  of  both  interference  sources;  therefore,  the  steepest 
descent  algor ithm  does  the  best  it  can  to  minimize  the  total  power. 

V.  CONCLUSIONS  AND  FUTURE  PLANS 

The  results  reported  here  demonstrate  that  a  useful  degree  of 
pattern  control  of  a  microstrip  adaptive  array  can  be  achieved  by  adjust¬ 
ment  of  the  reactive  loads  on  parasitic  elements.  The  specific  micro¬ 
strip  array  investigated  in  this  work  could  find  application  to  missile 
command  data  links  or  1 ine-of-sight  communications. 

The  main  advantage  of  a  reactively  steered  adaptive  array  (RESAA), 
as  compared  with  a  conventional  adaptive  array,  is  the  elimination  of 
the  mixers  and  other  hardware  needed  to  perform  the  complex  weighting 
of  the  output  of  each  element  at  an  intermediate  frequency  (IF).  In 
addition,  by  operating  without  the  need  for  an  IF  and  an  electronic- 
summing  element,  the  RF.SAA  can  probably  produce  a  larger  linear  dvnamic 
range  than  a  conventional  adaptive  array.  We  believe  that  a  RESAA 
configuration  can  be  made  with  substantially  fewer  components  than  a 
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FIGURE  24.  Beam  History  During  Steepest  Descent  Processing  for  Inter¬ 
ference  Incident  at  -90  Degrees.  For  this  history,  only  the  first  17 
iterations  of  the  50  needed  to  obtain  a  minimum  at  -90  degrees  are 
shown. 


TABLE  1.  Summary  of  Parameters  Used  in  Computations 
for  Figures  17  Through  25. 


Figure  Interference  angle 

number  of  incidence,  deg. 


17 

0 

0.01,0.035 

18 

-45 

0.01,0.028 

19 

-90 

0.01,0.025 

20 

0,-45 

0.01,0.28 

21 

0 

variable 

22 

0 

0.01 

23 

-45 

0.01 

24 

-90 

0.01 

25 

0,-45 

0.01 
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FIGURE  25.  Beam  History  During  Steepest  Descent  Processing  for  Two 
Equal  Amplitude  Interference  Signals  Incident  at  0  Degrees  and  -45 
Degrees . 


comparable  conventional  adaptive  array.  Another  advantage  is  the  neces¬ 
sarily  smaller  overall  size  (needed  to  maintain  tight  coupling  between 
elements) . 

There  are  two  disadvantages  to  the  RESAA  technique:  (1)  less 
control  of  the  pattern  is  obtainable  for  the  same  number  of  elements, 
since  only  one  control  device  is  used  for  each  element  with  the  RESAA 
compared  to  the  two  control  devices  (for  full  complex  weighting)  for 
each  element  of  a  conventional  adaptive  array  and  (2)  increased  control 
over  the  pattern  cannot  be  achieved  simply  by  adding  more  elements, 
since  new  elements  added  at  the  periphery  of  the  array  have  progressively 
less  influence  on  the  pattern. 

A  large  number  of  topics  relating  to  patlern  control  of  a  RESAA 
remain  to  be  investigated  and  these  studies  are  planned  for  l'V-81. 

These  topics  include  simulations  with  arrays  of  more  than  three  elements, 
nulling  of  more  than  one  interference  source,  the  effect  of  Gaussian 
noise  on  performance,  bandwidth  achievable  with  RESAAs,  the  equivalent 
number  of  degrees  of  freedom  that  can  be  realized,  and  the  practical 
aspects  of  implementing  an  experimental  RESAA  with  automatic  control. 
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Appendix  A 

measurement  OF  MUTUAL  impedance  matrix  elements  using 
AN  AUTOMATIC  NETWORK  ANALYZER 


Tiie  mutual  impedance  matrix  elements  required  for  Equation  2,  ,-t. 
are  difficult  to  calculate  for  microstrip  patch  antenna  elements. 
For  the  simulations  reported  here,  measurements  were  made  on  the  array 
shown  in  Figure  1,  and  these  values  were  used  in  the  simulation. 

A  Hewlett-Packard  (HP)  Model  8410  network  analyzer  controlled  by  an 
HP-9825  calculator  was  used  to  measure  the  scattering  parameters  of  the 
array.  Measurements  of  Sll  at  each  port  were  made  with  all  other  ele¬ 
ment  feed  ports  terminated  in  50-ohm  loads.  Measurements  of  S21  were 
made  pairwise  between  ports,  with  the  remaining  ports  terminated  in  50- 
ohm  loads.  These  measurements  of  Sll  and  S21  then  are  entries  in  the 
scattering  matrix  [ S ] . 

To  obtain  the  elements  of  the  impedance  matrix  [Z] ,  the  following 
equation  is  used:  *  3 

[Z|  =  [I  +  S)[I  -  S)"1  ,  ( A — 1 ) 

where  [I)  is  the  identity  matrix.  Table  A-l  lists  the  values  of  the 
elements  of  [Z]  obtained  in  this  manner  for  the  three  center  elements. 
These  values  were  used  in  the  simulations. 


TABLE  A-l.  Measured  Values  of  the  Mutual  Impedance  Matrix 
[ Z ]  Normalized  to  a  Resistance  of  1  ohm. 


Matrix  element 

Real  part 

Imaginary  part 

Z(11),Z(33) 

0.273 

-0.227 

Z(22) 

0.240 

-0.260 

Z(12),Z(21),Z(23),Z(32) 

0.138 

0.108 

Z  ( 1 3  )  ,  Z  ( 3 1 ) 

0.00780 

-0.0982 

13  R.  E.  Collin.  Field  Theory  of  Guided  Mn*es.  New  York,  McGraw- 
Hill,  1960. 
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Appendix  B 

RADIATION  I’AT'f KKNS  FOR  TOO-  AND  THREE - ELEMENT 
REACTI  V'EI.Y  STEERED  ADAPT]  VE  ARRAYS 


The  matrix  inverse  that  appears  in  Equation  8  in  general  cannot  be 
simplified  in  any  manner  that  shows  more  clearly  the  dependence  of  the 
radiation  pattern  on  the  reactive  loads.  However,  for  two-  and  three- 
element  arrays,  the  matrix  inverse  can  be  computed  explicitly.  This 
Appendix  gives  the  radiation  patterns  for  these  cases. 


TWO-ELEMENT  ARRAY 

For  a  two-element  array,  the  inverse  is  required  of  the  matrix: 


[ZA  +  V 


ZL1  Z12 
Z12  ZU  +  ZL 


(15-1) 


where  element-1  is  assumed  to  be  driven  and  element-2  has  a  reactive 
load  Ej  .  The  inverse  matrix  is  given  by 


|Z 


A 


+  Z,  1 


-1 


1_ 

D 


+  Z. 


(B-2 ) 


in  which  D  =  det(Z;\)  +  ZuZp-  For  an  element  separation  of  d,  the 
pattern  can  be  written  out  in  full  as 


F(0)  = 


i(°)  ’ 

D 


Z11  +  ZL  "  Z12e 


jk^dsinP-] 


where  g  (0)  is  the  element  pattern. 
3 


(B— 3 ) 
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THRF.K- ELEMENT  ARRAY 


The  inverse  is  required  of  the  matrix  given  by 


Z  +  L, 


IZA  +  V 


=  z. 


:  z., 


L 


Z1  Z2 


z  z. 

m  1 


Z1  Ze  +  4 


(B-4) 


where  Zt>  is  the  sel  f- impedance  of  the  end  elements,  ZIT,  is  the  self¬ 
impedance  of  the  middle  element,  Zj  is  the  mutual  impedance  between 
nearest  neighbors,  Z_>  is  the  mutual  impedance  between  the  end  elements, 
and  L'  and  L-j  arc  the  reactive  loads  on  element-1  and  -3,  respectively. 
The  inverse  is  given  hv 


Z  (Z 
m  e 

+V~Zl 

z1z2-z1zo-z1l3 

z.-z,z 

1  2  m 

|ZA  +  V1  -  B 

Z1Z2~ 

ZlVZlL3 

z2+z  L..+Z  L. 
e  e  3  el 

YW, 

T-LjEj-Z2 

'Yl 

zrz2 

z 

m 

ZlZ2-ZlZe-ZlL! 

Z  Z  +L.Z 
e  m  1  i 

(B-5) 


in  which  I)  is  the  determinant  of  the  matrix  [Z,  +  Z,  ] . 

,  A  L 

can  he  written  as 


The  pattern  then 


g  (0)  (  -jk  dsinO 

*  V  |Zm(Ze  +  S>  '  V* 


ik  dsinO 

+  Vz2  -  W  +  1Z1  -  w*  ° 


(B-6) 
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